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Abstract. The stability of trapped dilute Fermi gases against collapse towards 
large densities is studied. A hermitian effective contact-interaction for all partial 
waves is derived, which is particularly suited for a mean-field description of these 
systems. Including the s- and p-wave parts explicit stability conditions and 
critical particle numbers are given as function of the scattering lengths. The 
p-wave contribution determines the stability of a single-component gas and can 
substantially modify the stability of a two-component gas. Moreover it may give 
rise to a novel p-wave stabilized high-density phase. 
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Since the first realization of a Bose-Einstein condensate of ®^Rb atoms in 1995 ^ 
the field of trapped ultracold atomic gases experienced great experimental progress. 
This raised the question if a Fermi gas can be prepared under similar conditions and 
whether a transition to a superfluid state can be achieved. An important step towards 
a possible super-fluid state of trapped Fermi gases was the cooling of atoms to a 
temperature regime where degeneracy dominates 

A major problem for the evaporative cooling of fermions is the Pauli exclusion 
principle. The relative wave-function of two indistinguishable fermions has odd parity 
and hence they do not feel the s-wave part of the interaction which dominates the 
force between bosons at low kinetic energies. This limits the efficiency of evaporative 
cooling. Several techniques are under discussion to circumvent this problem, e.g., 
simultaneous trapping of two fermionic species j^, ^ , sympathetic cooling of a fermion- 
boson mixture [Q, and the use of p-wave resonances to enable efficient cooling via 
p-wave interactions ^, 0| . 

In this letter we investigate the question how p-wave interactions influence the 
properties of degenerate Fermi gases. The Thomas- Fermi approximation together with 
a new hermitian effective contact-interaction is employed to describe the effects of the 
atom-atom interaction in a dilute gas including s- and p-wave contributions. We use 
this formalism to investigate the influence of p-wave interactions on the stability of 
trapped one- and two-component Fermi gases against collapse towards high densities 
where the atoms escape from the trap due to three-body collisions and the formation 
of bound dimers. 

In the standard description of trapped atomic gases one uses point interactions 
in the s-wave channel only. In the following we present a hermitian effective contact- 
interaction (ECI) for all partial waves, which is derived to be an effective mean-field 
interaction. 

Consider two particles of mass m which interact via a spherical potential with a 
range much smaller than their relative wave length. An auxiliary boundary condition 
at a large radius leads to a discrete energy spectrum Eni , where n numbers the positive 
energy eigenstates and I denotes the relative angular momentum. Negative energy 
states, i.e., bound states, need not be taken into account since the ECI is used for the 
description of systems that are not self-bound. The energy shift AEni — Eni — Eni 
with respect to the eigenvalues Eni = k'^i/'m- of the kinetic energy without interaction 
can be expressed in terms of the phase shift r]i (k) for the Z-th partial wave § . 

For mean-field type calculations we require that the expectation value of the ECI 
calculated with mean-field states, i.e., free two-body states |nZrn/), equals the exact 
energy shift 

{nlmi \ \ nlmi) Ai?„;. (1) 

We choose the following hermitian ansatz for the operator of the ECI in a particular 
Z-channel 

The arrows above the derivatives indicate to which side they act. Using this ansatz 
in condition (|l|) we obtain an explicit expression for the interaction strength gi 
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where a; is the scattering length of the l-th partial wave. It turns out that (||) does 
not depend explicitly on the auxiliary boundary condition and hence we generalized 
the result to continuous momenta and energies. Approximation of the phase shifts 
r]i{k) in terms of the scattering lengths a/, so that gi is momentum independent, leads 
to deviations of at most 5% in the energy shift up to \kai\ ~ 3. For the s-wave and 
p-wave component this reduces to go = Anao/m and gi = 127raf/m, respectively. 

The s-wave part is identical to the widely used local contact interaction. Another 
way to include higher partial wave terms is the pseudo-potential of K. Huang and 
C.N. Yang 10|. This approach leads to non-hermitian interactions for I > 0. The 
two-body energy shift induced by this pseudo-potential is by a factor ^ii- smaller than 
the exact one, and is therefore not suited for mean-field calculations beyond s-waves. 

Using the s- and p-wave contributions of the ECI we calculate the Hartree-Fock 
energy-density functional of interacting Fermi gases composed of S distinguishable 
components in an external potential U{x). We use the Thomas- Fermi approximation, 
which is excellent for these systems [0 . The local energy density of the Fermi gas as a 
function of the local Fermi momenta = [Gn'^ p^{x)Y^^ of the different components 

^ = 1, . . . , S reads 

+ + E ^ 4 (^) + 4 (^)] ' (4) 

with coefficients 

^* ^ o?r^ — ' ^0 ^ fTT — ' -^1 ^ oTT"! — ■ ("^^ 

The first two terms are caused by the external potential and the kinetic energy, 
respectively. The remaining ones are the contributions of the s- and p-wave parts of 
the ECI. As a direct consequence of the Pauli principle, fermions of one kind interact 
only via the p-wave part, while fermions belonging to different components feel s- and 
p-wave interactions. The coefficient of the last term may be modified by effective range 
corrections of the s-wave channel. Since in the following applications oi is treated as a 
parameter this effect and others which contribute to the same order in are absorbed 
in fli. 

The ground state properties of the system are determined by minimization of 
the energy with respect to the local Fermi momentum under the constraint of a fixed 
particle number for each component. Implementing the constraints by means of 
chemical potentials fi^ the variational problem reduces to an algebraic equation for 
the local Fermi momentum. 

First we consider a single-component Fermi gas, where the extremum condition 
is given by 

nl[^i-U{x)]^l^^{x) + -^aU^ix). (6) 

Usually the s-wave contribution dominates at low energies, but since s-wave scattering 
is excluded here only p-wave scattering contributes to the mean-field energy. 

From the extremum condition (|^) we get the following upper bound for the local 
Fermi momentum of a single-component gas 

-ain{x)< , (7) 
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beyond which no energy minimum exists anymore. Condition can also be expressed 
as an upper bound for the chemical potential 

alm[f,~U{x)]<^^. (8) 

This limits the particle number of a stable single-component condensate for 
attractive p-wave interactions. 

The stability condition (|^) involves the local Fermi momentum or density at 
some point in space, which is not easy to access experimentally. Therefore, we 
express the stability condition by a phenomenological parameterization as function of 
particle number, scattering length and trap size. For simplicity we restrict ourselves 
to an external potential U{x) with the shape of a deformed harmonic oscillator, 
where £ = (lx£y£z)^^'^ = (m'^Wa^ajj^w^)"^/^ is the geometric mean of the oscillator 
lengths. The parameterization is fitted to the critical particle numbers obtained from 
the numerical solution of the extremum condition for a given scattering length 
using the maximum chemical potential (^) for a stable condensate. We find that the 
parameterized stability condition 

(v^y) - ^'^'^^^ ^ = -2.246, (9) 

for the single-component gas reproduces the numerically calculated stability limit 
within errors of 1%. 

The critical particle number iVc — [^/(Cai)]^ of a single-component gas for a 
typical set of experimental parameters (a\ = — 200aB, i ~ l//m) is about 8 x 10^. 
This is much larger than the particle numbers achieved in present experiments. 
Nevertheless, this stability limit could be exceeded in near future, e.g., by increased 
interaction strength or by the use of tightly confining traps . 

In the following we consider a gas composed of two fermionic species residing 
in two magnetic substates of hyperfine splitting. Although the s-wave interaction 
between the species usually dominates, we show that p-wave interactions cannot be 
neglected in general. 

The energy minimization using the energy-density functional for a two- 
component system leads to two coupled equations for ki [x) and ^2 {x) 

m[^ii - Ui{x)] = i fij{x) + ^0,0 kI{x) 

+ ^al[m4ix) + 3 4ix) + 5 4 ix)4 (x)] . (10) 

The second equation is generated by the exchange [^i — Ui{x)] ^ [/X2 — U2ix)] and 
Ki{x) <-> K2ix)- These two equations have a great variety of possible solutions, 
depending on the signs and the relative strength of Og and ai. In this letter we restrict 
ourselves to equal particle numbers or [fj, — U{x)] = [fii — Ui{x)] = [122 — U2{x)] which 
leads to equal local Fermi momenta k{x) — ki{x) — K2{x) and the single extremum 
condition for both components: 

12 4 

m[fi - U{x)] = H flo K^{x) H a? k^{x). (11) 

2 Stt Svr 

The generalization to different chemical potentials, different interaction strengths or 

more components is straightforward. 
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For the two-component system with s-wave interactions only, the stabihty 
condition was studied in |ll|. Using ( |ll| ) which includes the p-wave interaction we 
obtain the more general condition 

- aoK(f) - 2(aiK(x))^ < I (12) 

for each component not to collapse just because of mean-field effects. For ai — this 
reduces to the stability condition given in jl^, while for ai < the condition limits 
the particle number of a metastable two-component fermion condensate even more. 
It also shows, that a repulsive s-wave interaction (ao > 0) does not guarantee a stable 
condensate at larger densities if an attractive p-wave component is present. In this 
context it is also interesting to study the transition to spatially separated components. 

Several very interesting effects appear for interactions with attractive s-wave 
(ao < 0) and repulsive p-wave components (ai > 0). Here ai need not be the 
true p-wave scattering length but may include effective range and other effects which 
contribute to the K*-order in the energy density. In figure ^ the r.h.s. of the 
extremum condition (^) is depicted as function of the local Fermi momentum k. 
For an attractive pure s-wave interaction (lower curve) the r.h.s. exhibits an absolute 
maximum which leads to an upper bound for the chemical potential as discussed for 
the single-component gas. If a repulsive p-wave contribution is added, then the r.h.s. 
of dn] ) grows at large Fermi momenta due to the leading contribution. For values 
of the s- and p-wave scattering lengths which satisfy the condition 

^>:r47I~ 0-311 (13) 
|ao| Svr^/'^ 

the stability condition ( p^ ) is fulfilled for all values of the local Fermi momentum, i.e., 
the r.h.s. of the extremum condition is a monotonic function of k and the maximum 
does not appear anymore (upper curve). Thus the repulsive p-wave interaction 
stabilizes the condensate for any particle number even though the s-wave interaction 
is attractive. 

If the ratio of the scattering lengths is below the limit given by (|l^), then the 
r.h.s. of the extremum condition exhibits a low and a high density branch, as shown 
by the middle curve of figure [|. For chemical potentials fi below the value of the local 
maximum one obtains the regular low-density solution everywhere in the potential 
U{x). An example is shown in insert (a) of figure |l|. For larger fj, novel high-density 
solutions exist in areas of the trap where m[/Lt — U{x)] is larger than the value of the 
local maximum. As depicted in insert (b) the density shows a discontinuous jump from 
the outer low-density to the central high-density phase. This phase is stabilized only 
due to the presence of the repulsive p-wave interaction and occurs when the condition 
( p^ is violated. For values of /i in between the maximum and the minimum a Maxwell 
construction with equal pressure in the high- and low-density regime may be used to 
identify the equilibrium density. 

However, when ai/|ao| drops below 2/(37r'^/^) the density in the high-density 
regime grows rapidly so that the three-body recombination rate and thus the trap loss 
is not small anymore. If ai/|ao| < ■y/160/(7297r2) « 0.281 the local minimum of the 
r.h.s. of ( pr] ) occurs at negative values, so that even a self-bound solution exists. In 
this case the product |ao| k{x = 0) > |7r, which means that in the self-bound area in 
the center of the trap the mean distance p~^l^ between the atoms gets close to the 
scattering length |ao| and the underlying approximations break down. Nevertheless 
we expect that for interactions with ai/|ao| « 0.3 a metastable high-density phase 
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Figure 1. R.h.s. of equation ( |ll[ ) as function of the local Fermi momentum 
K. The tVire p cu rves correspond to different interaction parameters as labeled. 
Solutions of (hlh on the dotted parts of the curves correspond to maxima of the 
energy density, while the full branches denote (local) minima. The inserts (a) 
and (b) show the densities for two values of fj. indicated by the thin lines for a 
harmonic trap with mean oscillator length £. 



occurs, where the central densities are by one order of magnitude higher than in the 
outer low-density region. 

Like for the single-component gas we express the stability condition in terms of the 
particle number N of each component, the mean oscillator length ^ of a harmonic trap, 
and the scattering lengths oq and oi. We use the following slightly more complicated 
parameterization for the stability condition 

Co{^J)+Cf{^^)\c-,^'{^^){^^J < 1. (14) 

The fit has to be done separately for each type of attractive interaction. The resulting 
optimal coefficients are summarized in table ^ and reproduce the numerical stability 
limit again with deviations of less than 1%. 

The critical particle numbers A^c resulting from the stability condition (|lj) are 
shown in figure § as function of ao and ai. For negative p-wave scattering length Nc 
describes the maximum particle number of a metastable condensate. The white area 
in the figure marks the region where from the mean-field point of view a metastable 
condensate with smooth density exists for all particle numbers. In principle a positive 
ai always leads to a stable mean-field solution. But for oq < and values of 
'^i/l'^ol ^ 0.28 this solution is at a density which is too large for the metastable 
state. Therefore, the limit for iVc resulting from ( p^ is plotted. The bending over of 



Table 1. Parameters of the fitted stability condition ( |l4| ) for the two component 
condensate for different interaction types. 



interaction type 


Co 


Ci 


Coi 


n 


ao < 0, ai < 


-1.835 


-2.570 


0.656 


1 


ao > 0, ai < 


-1.378 


-2.570 


1.360 


1 


ao < 0, ai > 


-1.835 


-1.940 


2.246 
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the contour lines close to ai/|ao| ~ 0.3 indicates the onset of the novel high-density 
phase. 

As an example for estimating from figure ^ the critical particle number we take 
a two-component ^Li condensate in a trap with £. = lum. The extraordinary large 
attractive s-wave scattering length of oq — — 2160aB fl^, i.e. a^/i w —0.1, leads to 
the rather small critical particle number Nc — 26000 for each of the two components. 
As also seen from figure |^ the p-wave interaction will change significantly whenever 
|ai|/|ao| >0.2. 

We conclude that for Fermi gases p-wave scattering should not be neglected from 
the outset. If attractive it constrains the total particle number in the metastable 
condensate, if repulsive it helps to stabilize multi-component systems with attractive 
s-wave interactions. In special cases even a novel high-density phase may occur in the 
center of the trap which should not be confused with Bosc condensation of Cooper 
pairs. 
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